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Abstract 

Modern effective-theory techniques are applied to the nuclear many-body problem. A novel ap- 
proach is proposed for the renormalization of operators in a manner consistent with the construction 
of the effective potential. To test this approach, a one-dimensional, yet realistic, nucleon-nucleon 
potential is introduced. An effective potential is then constructed by tuning its parameters to 
reproduce the exact effective-range expansion and a variety of bare operators are renormalized in a 
fashion compatible with this construction. Predictions for the expectation values of these effective 
operators in the ground state reproduce the results of the exact theory with remarkable accuracy 
(at the 0.5% level). This represents a marked improvement over a widely practiced approach that 
uses effective interactions but retains bare operators. Further, it is shown that this improvement 
is more impressive as the operator becomes more sensitive to the short-range structure of the po- 
tential. We illustrate the main ideas of this work using the elastic form factor of the deuteron as 
an example. 

PACS numbers: 
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I. INTRODUCTION 



The construction of effective interactions for use in shell-model studies of nuclear struc- 
ture has enjoyed a resurgence due in part to the development of modern effective-field 
theories 0, || . Further, tremendous advances in raw computational power and numerical 
techniques have enabled the consistent and systematic implementation jiyji of 25-year-old 
approaches based on the so-called similarity-transformations methods |a lit la l| • Such im- 
plementations bypass most of the recent criticism levied on frequently employed shell-model 
approaches that rely on effective interactions that do not follow in any systematic way from 



a realistic nucleon-nucleon (NN) interactions [lfj, |TjJ. Indeed, similarity-transformations 
methods indicate how bare interactions and operators should be modified in a systematic 
way to account for the inevitable effects of truncations. 

Earlier work by two of us (JP and JRS) on low- dimensional quantum magnets [12, EH, 
3 2 3] made us familiar with a variety of theoretical approaches that have been recently 



adapted to the nuclear many-body problem (l7| . In that work it was shown how to combine 
similarity-transformation methods, specifically the Contractor REnormalization (CORE) 
approach of Morningstar and Weinstem|l8l. figj | . with effective interactions methods, such 
as those discussed by Lepage In particular, accurate predictions for the ground-state 
energy of the three-body system were made with relatively little computational effort when 
both techniques were used in a complementary fashion. As discussed in other recent publi- 
cations [i^, these similarity transformation methods may be understood in the context of 
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effective theories, which in turn rely on renormalization-group techniques 

What it is not at all clear (at least to us) in effective-theory approaches, is how to modify 
operators in a manner consistent to the modifications of the underlying Hamiltonian. The 
need for consistently modified operators must be emphasized. Parametrized operators are 
often added to improve quantitative agreement with data, but their origin is left unclear. 
Even in approaches based on similarity transformations where the modification to operators 
is well delineated, it remains common practice to employ bare (rather than renormalized) 
operators. In this work we adopt some modern concepts of low-energy Effective Theories 
(ET's) in the hope of improving some of these shortcomings. The basic assumption of ET's 
is that the complicated, and likely unknown, short-distance details of a theory are hidden 
from a long-wavelength probe. It should then be possible to modify the corresponding 
portion of the potential leaving its low energy properties intact. In order to achieve this, 
the low-energy properties must be known in advance either from experiment or, as in the 
case of this study, from solving the bare theory exactly at low energy. As has been observed 
by many authors 2^, 23{, ET's for the NN interaction reproduce the effective range theory 
of many decades ago. Part of the inspiration for the operator methods reported here arose 
from an especially simple derivation of effective-range theory which appears in the texts by 
Schiff and Taylor HQ. 

The paper has been organized as follows. In Sec. |H] a simple derivation of the effective- 
range expansion in one spatial dimension is presented. Next, a prescription for the renormal- 
ization of effective operators that is consistent with the construction of the effective potential 
is introduced. In Sec. IIHI expectation values for various effective operators are computed 
and are then compared to those obtained in the bare (exact) theory. Finally, conclusions 
and some ideas for the future are discussed in Sec. IIVI 
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II. FORMALISM 



The aim of this section is to adapt a textbook derivation of the effective-range formula 
in three spatial dimensions 24], HI to the one- dimensional problem considered here. These 
ideas are then used for the construction of an effective interaction that reproduces the 
scattering length and effective range of the exact (i.e., bare) theory. Finally, an approach 
is proposed for the renormalization of operators in a manner which is consistent with the 
construction of the effective interaction. 



A. Effective-range formula 

To arrive at the effective-range formula we proceed along the lines of Schiff and Taylor |2£ 
25j ] . adapting their derivation to the one- dimensional case considered here. The even-parity 
solution of the scattering problem satisfies the time-independent Schrodinger equation 



f- 2 +k 2 ~ U{x) 
dx z 



ip k (x) = (k 2 = 2fiE and U(x) = 2/iV(x] 



subject to the following boundary conditions: 
lim^ fc (x) = 1 + C(x 2 ) , 



x->0 



lim ip k (x) — 0& 0*0 = cos(kx) — tan. S(k) sin(A;x) 



cos (kx + 6(h)) 
cos S(k) 



(1) 

(2a) 
(2b) 



Note that <p k (x) denotes the solution of the free Schrodinger (U(x) = 0) equation that 
coincides with ip k (x) at large x. It then follows immediately from Schrodinger's equation 
that 



(I (■ 

dx 



k 2 %l) k (x)ip Q (x) 
k 2 (p k (x)(j)o(x) 



where the Wronskian of / and g is defined as 



fix) 






fix) 


g'(x) 





W(f,g)(x) = 

Upon integrating the difference of Eqs. (J3J) one obtains 

W(<f> k ,<f>o)(x)-W(il> k ,il> )(x 



f(x)g'(x) - f{x)g(x) 



(3a) 
(3b) 

(4) 



}r I dx[(j) k (x)(j)o(x) - xp k (x)i) Q (x) ) . (5) 



The contribution from the upper limit of the integral to the left-hand side of the equation 
vanishes, as ip k (x) =4> k (x) at large distances. Further, as the derivative of the exact scattering 
solution vanishes at x = [see Eq. (|2aj)] the Wronskian W(ip k , ipo) vanishes as well. This 
yields 



TU(0 fc ,0 o )(x = O) = fc (O)^(O)-^(O)0o(O) 



k l I dx((f> k (x)(f) (x)-i/; k (x)il) (x)) , (6) 
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which in turn generates the well known effective-range formula 



kttm5(k) = — — k 2 [ dxU k {x)^ {x) - Mx)M%)) = — - T ^k 2 + O(k^) . (7) 
ao Jo v / a 2 

Note that the (even-parity) scattering length and effective range parameters have been 
defined as 



a, 



oo 

1 - limfctan<J(A;) , r = 2 / dx(ip 2 (x) - <j>l(x) ) . (8) 
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B. Effective Potential 

The purpose of this section is to summarize briefly the main points from Ref. 0] which 
will, in turn, motivate our proposed method for constructing effective operators. To start, a 
bare one-dimensional NN interaction with the same patholog 1GS £lS du realistic interaction is 
assumed. That is, the bare potential is given by the sum of a strong short-range repulsive 
and a medium-range attractive exponentials: 

V(x) = V s e- m °W + v y e" mv|a;| . (9) 

The two masses were chosen to be equal to m s = 400 MeV and m v = 783 MeV, respectively, 
while the strengths of the potentials (V s = —506 MeV and V v = +1142.49 MeV) were cho- 
sen to give a binding energy and point root-mean-square (rms) radius for the symmetric 
("deuteron") state of E h = — 2.2245 MeV and r rms = 1.875 fm, respectively. 

Employing an option originally suggested by Lepage [l|, and later adapted by Steele and 
Furnstahl jia H?} to treat the iViV interaction, we propose a gaussian cutoff for the effective 
potential of the form 

1 ( d 2 \ 
V eS {x) = -(c + d— + e— + ...\exp{-e)-, i = x/a. (10) 

The parameters of the effective potential (c, d, e, . . .) are fixed to reproduce the low-energy 
scattering phase shifts. That is, one adjusts the parameters until the following equation is 
satisfied: 



A;tan<5(&;) = k 2 [ dx(ipk{x)ipo{x) — 4>k{x) 4>o{%)) (H) 

ao Jo V / 

= 1 _ k 2 [°° dxUf(x)if>f(x) - faWMz)) , (12) 

dn In V / 



where ^l s (x) is a scattering solution of Eq. (JTJ) with V(x) replaced with V e s(x). Note that 
as in Ref. [l7J, the gaussian cutoff parameter has been fixed at a = 1.16 fm. The above 
condition may be rewritten in the following convenient form: 

5(T)(k;c,d,...) = J dx(il) k (x)ip Q (x)-i>f(x)i>f(x)) = 0. (13) 

Evidently, it is neither demanded nor expected that Eq. (|T3~j) be satisfied for arbitrary large 
values of k. Rather, one follows a hierarchical scheme, based on power counting, that assures 
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that observables in the bare and effective theory be indistinguishable at low energies. It 
should be emphasized that the specific form of the potential is somewhat arbitrary, as the 
short-range structure of the theory becomes encoded in the effective parameters. 

As a simple illustration of this procedure we display in Fig. ^ bare (solid line) and ef- 
fective (dashed line) NN potentials. The bare potential, with its characteristic short-range 
structure, yields a scattering length of a = 5.247 fm and an effective range of r = 1.521 fm, 
respectively. The calculation of low-energy phase shifts is repeated using the effective poten- 
tial [Eq. (fTUj) ] with its two parameters (c and d) adjusted to reproduce the exact effective- 
range expansion to order k 2 . The resulting parameters (c = —0.039 and d = —0.160) are 
natural and yield a smooth effective potential which as far as the low-energy properties of 
the theory are concerned, is practically indistinguishable from the bare potential. Indeed, 
bulk properties of the ground-state (henceforth referred to as the "deuteron" ) are predicted 
to be identical to those obtained in the bare theory. This in spite of the vastly different 
short-range structure of the wavefunctions (see inset in Fig. Q). 
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FIG. 1: Bare (solid line) and effective (dashed line) NN potentials in a realistic 1-dimensional 
model. The inset shows "deuteron" ground-state wavefunctions. The sharp features of the bare 
potential are no longer present in the effective potential. Although the short-distance structure of 
the wavefunctions are different, the exponential falloff (binding energy) is unchanged. 



C. Effective Operators 



One of the main criticisms levied on traditional shell-model approaches is the lack of 
consistency between the construction of the effective potential and the renormalization (if 
any) of the bare operators 0, 11 1. While important steps have been taken to correct this 
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inconsistency, both in the area of low-dimensional quantum magnets [14| and nuclear struc- 
ture Q , these are in the very early stages. Further, it is unknown how to construct effective 
operators that are consistent with both the effective theory and similarity-transformation 
based approaches. 

In this contribution we propose an approach for the modification of operators that is 
consistent, indeed mimics, the construction of the effective potential. We assume that the 
momentum dependence of matrix elements of (simple) single-particle operators may be ac- 
counted for by an expansion in powers of k 2 having the same form as the effective-range 
formula [Eq. (JJJ)]. To do so, one demands that matrix elements of effective operators (C e ff) 
with scattering-wave solutions of the effective potential ("01 ) possess the same momentum 
dependence as those using the bare operators with the exact wavef unctions. In analogy with 
the definition of the effective potential [Eq. (|10|)]. we parametrize the effective operators via 



O eS {x;c,d,...) = 0{x) 



,2 



1 + ( c + d— - + ...) exp(-£ 



(14) 



This parameterization affects only the short-range behavior of the operator just as using 
the effective potential modifies only the short-range structure of the wavefunction. The 
parameters c, d, . . . (as before) are tuned to the low-k 2 behavior of the exact matrix elements. 
To be more specific about our procedure, we fit — in complete analogy to Eq. (JT3j) — the 
parameters of the effective operator by requiring that 

6(0)(k;c,d,...) = J dx^M^)0(x)^(x)-^f(x)O eS (x;c } d 7 ...)^(x)^ = 0. (15) 

The integral in this expression is convergent as the effective theory demands that 

lim ijjf (x) = ipk(x) , (16a) 

x— >oo 

lim O eS {x; c, d, . . .) = 0{x) . (16b) 

However, to insure that each term separately is convergent, we add and subtract the following 
term: 

dx{4> k {x)0{x)^{xYj , (17) 



where we recall that (j) k {x) is the free solution of the ID scattering problem [see Eq. (J2| 

To extract the parameters of the effective operator we now fit — in the spirit of the 
effective-range expansion — the low-energy matrix elements of the bare operator between 
bare scattering wavefunctions according to 

(O)bb(AO = f o dx(^ k {x)0{x)^{x) - fc (x)O(x)0^(x)) = a + (3k 2 + . . . (18) 

The parameters fixing the effective operators are then adjusted so that the above expansion 
is recovered. That is, 

oo 



(0) EE (k) = dx(if>f (x)O eS (x;b,c,...)^(x)-(/ )k ^)0(x) ( f) k ^(x)j =a+(3k 2 +... (19) 
Note that when 0(x) = O c g = l one recovers the effective-range expansion. 
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III. RESULTS 



In this section we compute matrix elements of various operators using three different 
schemes. The first scheme uses bare operators with bare wavefunctions (we label these 
calculations as "B+B"); these should be regarded as "exact" answers. Second, we compute 
matrix elements in an approximation (labeled as "B+E") that uses effective wavefunctions 
but retains bare operators. As we show below, for operators insensitive to short-range 
physics this inconsistency introduces small discrepancies. However, the more important the 
short-range physics, the greater the lack of accord. Finally, we perform calculations in a 
consistent low-energy approximation ("E+E") that employs both effective wavefunctions 
and effective operators. Showing that these calculations are in excellent agreement with the 
exact (B+B) answers represents the central result of the present work. 

Because of its simplicity, a natural place to start testing the proposed approach is the 
calculation of the root-mean-square radius of the deuteron, which is given by 

(x 2 )bb = (lj2 ( x " - x <™) 2 ) = r dx T^U x ) = (1-87977) 2 fm 2 . (20) 

Z n=l J -°° 

The corresponding calculation in the effective theory requires a renormalization of the bare 
operator. To do so we follow the prescription outlined in the preceding section [see Eqs. (JTHJ) 
and (fTTJj)] to obtain 



X cS — x 



d 2 

c + d—— ) exp(-£ 2 ^ 



. )r2 j ...... . , .. 1.520, d= -0.305) . (21) 

In this manner the root-mean-square radius predicted by the effective theory becomes 

f°° r 2 / \ 2 
(x 2 ) EE = dx^fU^(x)) = (1.87988) 2 fm 2 . (22) 

J —oo ^ 

This represents a discrepancy of about 1 part in 10 4 . While this result is gratifying and lends 
some credibiliity to the approach, it hardly qualifies as a stringent test of the formalism. 
Although both the effective operator and the ground-state wavefunction are modified at 
short distances (see Fig. Q and 121) the operator itself has so little support at short distances 
that the two integrands [Eqs. (J20j) and (J22)) ] become practically indistinguishable from each 
other (see inset on Fig. |2J. Indeed, an acceptable result is obtained even when the operator 
is not properly renormalized: 

f°° T 2 / \ 2 

(x 2 ) BE = dx-U^(x)) = (1.87834) 2 fm 2 . (23) 

J — oo ^ 

A more sensitive test of the approach is provided by the elastic form factor of the deuteron, 
which in our simple one dimensional model reduces to the following expression: 

F el (q) = W)? , P{q) = r dxcoB (?) = 1 - tr{x 2 ) + 0(q 4 ) . (24) 



2 J rgs ^ y " 2 

The corresponding expressions in the B+E and E+E approximations are given by 

p(g) BE = dx cos (T) (4>£{xj) 2 , (25a) 



-oo 
oo 



2 



(>Ui)ee= t <l-r <"^(f)j eff • (*)1>) 
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x(fm) 



FIG. 2: Bare (solid line) and effective (dashed line) X 2 operator. Note that while the operators are 
considerably different at short distances, its ground-state expectation values are not (see inset). 



with the effective operator renormalized at short-distances as detailed above. That is, 



cos 



qx 

Y 



cos 



J eflf 



qx 

y 



1 + ( c{q) + d{q)— j oxp( 



(26) 



The renormalization procedure is illustrated in Fig. El at the single momentum-transfer value 
of q = 2 fm~ . It is important to note that the renormalization coefficients (c and d) must 
be tuned at each value of the momentum transfer q. The solid line in the figure displays 
an effective-range-like expansion for the bare operator 0(x) = cos(qx/2) [as described in 
Eq. (|18|) . with the slope and intercept clearly indicated in the figure (note that for clar- 
ity the plots are normalized to one at k 2 = 0). It becomes immediately evident that the 
predicted low-energy behavior of the exact theory can not be reproduced without a proper 
renormalization of the operator. Indeed, the B+E calculation predicts the wrong momentum 
dependence; the sign of the slope is wrong! In contrast, it becomes a simple matter to tune 
the parameters of the effective operator to reproduce exactly the low-energy behavior of the 
exact theory (dashed line). Having corrected the short- distance structure of the operator 
one proceeds to compute the elastic form factor of the deuteron, which now is a prediction 
of the effective theory. The structure of the form factor (again at q = 2 fm _1 ) is shown in 
Fig. |U The main panel shows bare (i.e., cos(x)) and effective operators, displayed as solid 
and dashed lines, respectively. Both the effective deuteron wavefunction (inset on Fig. ^) and 
the effective operator differ considerably from their bare counterparts at short distances — 
and so is the product of the (square of the) wavefunction times the operator (inset on Fig.0J). 
Yet the area under the curve — whose square is proportional to the elastic form factor — is 
essentially unchanged. For comparison, the exact (B+B) and effective (E+E) theories yield 
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FIG. 3: An effective-range-like expansion for the elastic form factor of the deuteron at q = 2 fm~ . 
Calculations are displayed for the bare theory ("B+B"), the effective theory ("E+E"), and for a 
"hybrid" approximation that uses bare operators with effective wavefunctions ("B+E"). The solid 
(B+B) and dashed lines (E+E) are identical (by construction) since the effective parameters (c 
and d) are tuned to reproduce the expansion for the bare theory. 



values of F el = 0.02019 and F e i = 0.02012, respectively. In contrast, a (B+E) calculation with 
an effective wavefunction — but still employing a bare operator — results in a discrepancy of 
nearly 20% (F e \ = 0.017062). We conclude the discussion of the elastic form factor of the 
deuteron by displaying in Fig. EJits momentum-transfer dependence up to q = 5 fm -1 . Recall 
that effective parameters must be tuned for every value of q. It is evident from the figure 
that the renormalization of the operator at high-momentum transfers is essential, as it is at 
high q that the short- distance structure of the wavefunction (and of the potential) is being 
probed. Failing to correct the operator results in a rather poor representation of the elastic 
form factor for q > 2 fm -1 (squares). 

For completeness, and as a further stringent test of the formalism, we compute ground- 
state observables for an operator with the most extreme short-range structure possible: the 
Dirac delta function. In the bare theory the ground-state expectation value is simply given 
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x(fm) 

FIG. 4: Bare (solid line) and effective (dashed line) for the cos(x) operator. The inset shows the 
product of the square of the wavefunction with the operator for the three calculations discussed in 
the text. The elastic form factor of the deuteron (at q = 2 fm^ 1 ) is proportional to the square of 
the area under the curve. 



by the square of the deuteron wavefunction at the origin. That is, 

/oo 
dx6(x)ip^(x) = V&(0) = 0.01521fm~ 1 . (27) 
-oo 

As the sharp features of the bare potential are no longer present in the effective potential, 
the effective deuteron wavefunction at short distances is considerably larger than the bare 
wavefunction (see Fig. GJ. As a result, a (B+E) calculation using an effective wavefunction 
but a bare delta-function operator grossly overestimates the result: (5(x))BE = 0.04981fm _1 . 
Instead, by following the renormalization procedure outlined above, one obtains an effective 
Dirac delta-function operator, 

1 / $2 \ 

5 eS {x) = - f c + d-^j exp(-£ 2 ) (c = 0.250, d= -0.035) , (28) 
that yields a ground-state expectation value of (S(x))ee — 

0.01523 fin -1 . This result deviates 

from the exact value by less than one part in a thousand. 

In Table IJ we have listed (for completeness) some of the results presented previously in 
graphical form. The operators appearing in this table are listed in order of the importance 
of their short-range components. For example, the root-mean-square radius of the deuteron 
depends little on the short-range structure of the wavefunction while the Dirac 5-function 
operator depends exclusively on it. For each operator, we show the two effective coefficients 
c and d determined by the fitting procedure outlined above. All these are dimensionless 
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FIG. 5: The elastic form factor of the deuteron in the three calculations discussed in the text. Note 
that the predictions from the effective theory (E+E) agree with the exact theory at at all values 
of the momentum transfer q. 



quantities and it is gratifying that they are all of order one in keeping with the priciple of 
"naturalness" [1^]. We note that in some cases failing to renormalize the operator (B+E) 
leads to discrepancies that are as large as 50% or even 100%. In contrast, calculations using 
effective wavefunctions and effective operators (E+E) show excellent agreement with B+B 
calculations regardless of the short-range structure of the operator. We stress that all these 
are predictions of the effective theory, as the tuning of parameters is done in the scattering 
sector. In particular, it is satisfying that the elastic form factor of the deuteron at q = 
deviates from unity by less than 2 parts in a thousand. We emphasize that such precise 
agreement is non-trivial; indeed, it reflects the soundness of our method. 



IV. CONCLUSIONS 



While the field of nuclear structure has benefited from recent advances in numerical 
algorithms and sheer computational power, the shell-model problem, in its purest form, 
remains intractable. As a result, an important part of the nuclear-structure program for 
many years has focused on the construction of effective interactions for use in shell-model 
calculations. Two of the most promising approaches are based on the so-called similarity 
transformation methods (in its many varieties) and on effective-field-theory techniques. The 
main tenet underlying both approaches is that the short- distance structure of a theory 
(which is complicated and at present unknown) are hidden to a long-wavelength probe. It 
should then be possible to "soften" the corresponding short-range portion of the potential 
while leaving all low-energy properties of the theory intact, thereby providing a significantly 
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(0(x)) 


c 


d 


B+B 


B+E 


E+E 


V< x 2 > 


+1.51979 


-0.30545 


1.87997 


1.87834 (0.09%) 


1.87988 (0.01%) 


|(cos(0))| 2 


-0.02016 


+0.02985 


1.00000 


1.00000 (0.00%) 


1.00155 (0.16%) 


|(cos(x/2))| 2 


-0.06440 


+0.13429 


0.04159 


0.04261 (2.45%) 


0.04181 (0.53%) 


|(cos(x))| 2 


-0.79386 


-0.54503 


0.02019 


0.01706 (18.4%) 


0.02012 (0.35%) 


|(cos(2x))| 2 


-0.52925 


+0.14806 


0.01541 


0.00698 (220%) 


0.01544 (0.19%) 


(5(x)) 


+0.24964 


-0.03501 


0.01521 


0.04981 (327%) 


0.01523 (0.13%) 



TABLE I: Ground-state expectation values for various operators in the different approximation 
schemes described in the text. B+B indicates bare operators with bare wavefunctions, B+E bare 
operators with effective wavefunctions, and E+E effective operators with effective wavefunctions. 
Note that the root-mean-square radii are given in fm and (5(x)) in fm _1 . Finally, the c and d 
coefficients are dimensionless parameters of the effective theory. 

more tractible - from a computational point of view - interaction. 

The main focus of the present paper is the determination of single-particle operators 
which can be employed consistently in conjunction with wavefunctions obtained using ef- 
fective interactions. As observed by many authors, such consistency is essential to the 
correct implementation of effective theories. For computational simplicity we adopted a 
one-dimensional NN interaction, that nevertheless incorporates the well-known pathologies 
of a realistic AW potential. The central result of this work is the proposal and implemen- 
tation of a single underlying approach for the construction of both effective interactions 
and effective operators. The construction of the effective interaction follows a well-known 
approach that is based on a textbook derivation of the effective-range expansion. What is 
not well known (at least to us) is that the same approach may be generalized to effective 
operators. 

Results from such an implementation are very gratifying, as evinced from a variety of 
calculations of ground-state observables. For those observables insensitive to the short-range 
structure of the potential, such as the root-mean-square radius of the deuteron, the renormal- 
ization of the bare operator, while required by consistency, is of little numerical consequence. 
Yet, failing to properly renormalized operators sensitive to short-range physics, such as the 
elastic form factor of the deuteron at high-momentum transfers, can yield discrepancies as 
large as 200%. The consistent renormalization procedure advocated here yields in all cases 
errors of less than 1%. 

We conclude with a short comment on future work. The results presented here are 
encouraging and lend validity to the proposed approach, which is currently being extended 
to the three-body system. Different algorithms are being employed to solve for the ground- 
state of the three-body system and in all cases, perhaps not surprisingly, better convergence 
properties are obtained with the effective rather than with the bare interaction. The results 
obtained here also constitute a promising first step toward our ultimate goal of combining 
similarity-transformation methods with effective interactions. The effective interactions and 
operators obtained here — with their sharp short-range features no longer present — could 
provide a more suitable starting point for the numerically-intensive approaches based on 
similarity transformations. Finally, the method proposed here will have to be extended to 
three-spatial dimensions. Other than numerical complexity, we do not foresee other serious 
challenges. Indeed, the approach presented here for the construction of effective interactions 
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(whose three-dimensional derivation may be found in several textbooks) had to be adapted 
to one spatial dimension. In summary, a novel approach for the renormalization of operators 
in a manner consistent with the construction of the effective potential has been proposed 
and implemented with considerable success. The results obtained here are gratifying and 
suggest how in the future effective theories may be profitably combined with more traditional 
methods to tackle the nuclear many-body problem. 
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